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J> ' Abstract 

In this paper we investigate supersymmetric effects to the threshold production 
cross section of top quark pairs in electron positron annihilation. In particular, we 
consider the complete one-loop corrections from the strong and weak sector of the 
. Minimal Supersymmetric Standard Model. 
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1 Introduction 



One of the main goals of a future electron positron collider is the precise measurement of 
the top quark production cross section in the threshold region. The comparison to the 
theoretical prediction allows for a precise extraction of the top quark mass, its width, the 
strong coupling and — in case the Higgs boson is not too heavy — the top quark Yukawa 
coupling. 

The next-to-next-to-leading order (NNLO) QCD corrections to the total cross section 
a(e + e~ — > tt) has been completed several years ago [1]. One observes large perturbative 
corrections from the second order terms which make a precise prediction difficult. In the 
recent years a big effort has been undertaken to complete the third-order corrections to 
<r(e + e~ — > tt) [2-12] . First numerical estimates [13] indicate that the convergence of the 
perturbation theory is improved after the inclusion of the NNNLO terms. In addition to 
the third-order corrections also the resummation of the NNLL terms is studied [14-16]. 

In order to profit from precise experimental measurements it is desired to reach an 
uncertainty below approximately 3% from the theory side [17]. Radiative corrections 



of this order can easily be reached by effects from theories beyond the Standard Model 
(SM). In this paper we consider the effect of supersymmetric corrections within the Mini- 
mal Supersymmetric Standard Model (MSSM). Furthermore we confirm the results from 
Refs. [18-20] and [21] obtained in the framework of the SM and two-Higgs-doublet model 
(THDM) of type II, respectively. 

It is convenient to perform the calculation of the production cross section in the 
framework of an effective theory where the produced top quarks are described by a non- 
relativistic two-particle Green's function. All effects connected to energy scales above 
fi ~ mw are contained in coefficient functions which represent the new couplings in the 
effective Lagrangian. Since the masses of the supersymmetric particles are above the 
electroweak scale they only influence the matching coefficients of the effective operators. 
For the top quark production we have to consider the vector current in the full and 
effective theory which constitutes a building block for all threshold phenomena involving 
the coupling of the initial electron and positron via photon, Z boson or box diagrams to 
heavy quarks. 

The remainder of the paper is organized as follows: In the next Section we provide the 
formulae which are necessary for the evaluation of the threshold cross section. Afterwards 
we discuss in Sections [3] and H] the numerical effects from the strong and weak sector of 
the MSSM and present in Section 5 our conclusions. 

2 Framework 

Non-relativistic QCD (NRQCD) allows for a consistent separation of the hard corrections 
connected to energy scales of the order of the weak gauge bosons or higher from the 
soft scales which are involved in the top anti-top boundstate. Within NRQCD we can 
normalize the production cross section to cr(e + e~ — > = (47ra 2 )/(3s) and denote 

the ratio by R 

R(e + L e~ R ^ttX) = ^-lm[(h Ry ) 2 H v + (h RA ) 2 H A }, (1) 

where s is the square of the center-of-mass energy. In Eq. (TjQ) left-handed positrons and 
right-handed electrons are considered; for e^e^ in the initial state a similar expression 
is obtained by replacing R by L in Eq. ([]]). Note that the initial states e^e^ and e\e^ 
are suppressed by a factor {m e /M w ) 2 ~ 1CT 10 and are thus negligible. h R y an d h R A 
are so-called helicity amplitudes which absorb the matching coefficients representing the 
coupling of the effective operators. They take care of the hard part of the reaction. The 
first subscript of h refers to helicity of the electron, and the second one to the vector 
(Jy = ■07 /i ^) or axial-vector coupling (J^ = rfj^jsi/j) of the gauge bosons to the top 
quark current. In this paper we evaluate corrections to h RV and h L V . 

The bound-state dynamics is contained in the so-called hadronic part formed by 
current-current correlators within NRQCD. They are denoted by Hy and Ha in Eq. ([1]) 
and will not be considered further in this paper. At threshold the contribution from 
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the axial-vector current is suppressed by two powers of top quark velocity thus we only 
consider the vector current J v in this work. Its counterpart in the effective theory reads 

It is convenient to separate the photon and Z contribution in h^y and write 

hf$ = h]^ ee + h^ ee . (2) 

Here the tree-level contributions are given by (I = L/R) 

t^I V = QeQt j 
^Z,tree s Pj Pv 



s-MY 

t i at 



Pv 

P{ = W/z-W/ ^ (3) 



2 

(T 3 ) fl -s 2 w Q f 



where the pj is the coupling of a fermion (/ = e, t) to the Z boson, s w is the sine of the 
weak mixing angle (c^ = 1 — = m^,/m|), and electric and iso-spin charges for top 
quark and electron are given by 

Q e = -1, Q t = 2/3, (T 3 ) tL = l/2, (T 3 ) ei = -l/2, (T 3 ) fR = 0. (4) 

In the following the abbreviation T( = (T 3 ) f L will be used. Let us note that h^A can be 
obtained by substituting (3 V by P\ = {P l R — P t L )/2 in formula ([3]). The loop corrections 
are taken into account via 

hi, v = hTj + hfy, (5) 



where X stands for QCD, SQCD (supersymmetric QCD), SM, THDMll or MSSM. The 
numerical effects are discussed for the quantity 

A x = S& = 2^-Re (hl v ) + 2^-Re (h* v ) 
RL ° (hfy) 2 + (hf v ) 2 

where the sum over all helicity states of the incoming electron and positron has been 
performed. Let us note that in our case for the evaluation of hiy one has to set s = 4mf . 
Furthermore the external top quarks are on their mass shell. In addition we are only 
interested in hard corrections resulting from the real matching condition. Corrections to 
the cross section stemming from imaginary part of the matching coefficient, which takes 
into account the finite lifetime of the top quark, are discussed for SM in Ref. [22]. 



1 In this paper we use the THDM type II where u/d-type quarks couple to different Higgs doublets 
H u /Hd- Note that the Higgs sector of the MSSM corresponds to the Higgs sector of THDM type II. 
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Figure 1: QCD and SQCD diagrams, (a): Gluon contribution at the tt7/Z-vertex. (b): 
Gluon contribution to the top-quark selfenergy. (c): Gluino contribution at the tij/Z- 
vertex. (d): Gluino contribution to the top quark selfenergy. 



For the generation of the Feynman diagrams we use the Mathematica program 
FeynArts [23]. The amplitudes are further processed with the help of the programs 
FormCalc [24] and FeynCalc [25] which take the traces, map the occuring integrals to a 
standard basis and reduce the tensor integrals to a minimal set of scalar integrals usually 
denoted by Aq, Bq and Co- Since we have a quite particular momentum configuration it is 
not possible to use the above mentioned packages as black boxes but apply some modifi- 
cations. In fact, the choice s = 4m 2 allows for a partial fractioning in the denominators of 
the loop integrands appearing in tt^/Z- vertex and box diagrams which effectively reduces 
the number of external legs by one. Consider, e.g., the integrand of a generic three-point 
function (omitting the ie prescription) 

(7) 



2 + 2q x p - Mf + m 2 )0 2 - 2q 2 p - Mf + m 2 )(p 2 



Mf) 



where p is the integration momentum and q\ = q\ = m 2 are the squared momenta of the 



top quarks. After choosing q\ 
at 



q 2 = q/2 and applying a partial fractioning one arrives 



2 



Mf - m 2 



Mf+Mf . 2 

1 2 2 + mf 



Ml 



X 



+ 



p 2 + q ■ p — M\ + m 2 p* — q ■ p — iv± 2 



M 2 + m? t 



As a consequence the result can be expressed in terms of only two-point functions. In a 
similar way one can express the box diagrams in terms of three-point functions. 



3 Supersymmetric QCD 

In this Section we consider the effects from Supersymmetric QCD (SQCD) to the top 
quark threshold production. There are only four contributing Feynman diagrams which 
are shown in Fig. HJ The one-loop QCD corrections are known since long [26] and the 
corresponding matching coefficient is defined via the relation 
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This effect can be incorporated in the helicity amplitude by a simple rescaling of the 
tree-level contributions 



h? C v D = h^a-y + hf^a*. (10) 

Using the explicit expressions for h^y, one can already see that the contribution of a\ 
to the relative correction of the cross section is in general suppressed by factor 0.08 
compared to the one resulting from contribution of a? x . For the QCD the coefficients a?J Z 
read 

a"J Z = <£>. (11) 

Since we work in a supersymmetric framework we repeated the calculation of dp within 
dimensional reduction [27]. Although both the one- loop vertex corrections and the result 
for the wave function counterterm are different from their counter parts in dimensional 
regularization we observe that the result given in Eq. (Q does not change. This is expected 
since at tree-level the strong coupling constant is absent. 

The SQCD corrections can also be cast in the form of Eq. ( fTUl) 



^sqcd = f qp* < q + h z,*» a z ) (12) 

with 

4 = r^ + ^j-g^rg^, (13) 

where Ty~ represents the gluino contribution to the vector part of the one-loop vertex 
normalized by the corresponding tree level coupling. The wave function renormalization 
constant Z v j A - — 1 + 5Z V , A ~ is defined in the on-shell scheme and renders the jti and 
Zti vertex finite. The definition of both counter terms can be found in Appendix lAl 
where they are expressed in terms of vector- (V) and axial-vector part (A) of the top 
quark selfenergy. The subscript g reminds that only the diagrams involving a gluino 
are considered in each expression. Since the results are quite compact we present the 
analytical formulae for the individual contributions of the right-hand side of Eq. (fl"3l) . 
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The contributions to the wave function counterterm reads 



o.s 



s=l 



6nmi 



2m 2 



2m g m t (VL sls2 + fi S 2si) 



+ K + 



m-„ — m 



+ (Q slsl + Q s2 s2) A (m~) - A Q (mf) 



+ ( 



2 2 2 

m~ — m, — m~ 



) ^slal + ^s2s2^B ( 



2 2 2 
h, m ~g, m t s J I ' 



s=l * ^ 



- { m t + m l~ m l) B o W, m 2 g, ml) \ , 



and the vertex corrections are given by 



r 7 

1 V,g 



E 



a: 



:(^slsl + ^s2s2) 



2m 2 (m 2 - mfj B (4m 2 



2 2 



V,g 



^ 9nm 2 (m 2 + m| — m| ) 
+ | (m 2 + mf - m|) A (mf) - A (m|) + 2m 

~ \ \ m t ~ 2 ( m ? + m D m i + ( m l ~ m tY} B o m i m t) 

2 

E 



x 



+ 



+ 



ulsl ~r i'u2s2i 

i 97rm 2 {m 2 + m~ — 2m? — 2m 2 j (8^^ — 3) 

4 o / 2 i 2\ 2 i / 2 2\ 2 1 T) I 2 2 2 

i s ~ 2 \ m t + m g) m t s + \ m g ~ m t) B o {m t , m- g , m- t 



+ 



X 



+ m l - 2 ( m t + m l) 

A (m 2 ) - A (ml) + 2m 2 t (l + B (4m 2 , m|, m?J) 
%\ (ml +m 2 iu +2 (m 2 - m 2 )) - \ (m| - m?J 

Bo (m 2 t , ml, \m\ s + \m\ - m 2 ) j + j, 



(14) 



(15) 



In Eqs. ( TT4l) and ( fT5l) we introduced the abbreviations Qijki = U^U^ where Uij are 
the elements of the top squark mixing matrix (cf. Appendix [B]). The conventions for the 
functions Aq and Bq are adapted from Ref. [24, 28] where explicit results can be found. 
Further B' is the defined as derivative of B with respect to the first argument. Our 
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analytic formulae are in agreement with Ref. [29] where the result has been expressed in 
terms of a one-dimensional integral assuming a real mixing matrix for the top squarks. 

It is instructive to consider the limit where all SUSY particles have a common mass 
^susy- In this limit the above formulae are simplified significantly. In particular, the 
result becomes independent of the matrix elements Z7y and 5Z f A ~ = 0. We furthermore 
have a g = a~ = a~ which reads 

/ 2 \ J 1 2 3 4 d/ / 2 2 2 \ 

a sl m SUSYJ — Qm^TC 1 2 m * ol m t' m SUSY! m SUSYj 

+ ( m susv - m t ) [Bo( m t, ^susy, ™susy) - B (Am 2 t , m^ USY , m^ USY )] J 

= ^{y^W + W + M + o^)}- as) 

After the second equal sign we have expanded the result in terms of y = ^j/^susy- 

Let us in the following discuss the numerical effects of the one-loop QCD and SQCD 
corrections. For m t = 173.1 GeV and or s (mt) = 0.108 the QCD corrections amount 
to A QCD = 18.3% (corresponding to af\mz) = 0.1176) and thus constitute the largest 
contribution. 

In the simplified scenario described by Eq. ( 1T61) one obtains the SQCD corrections as 
shown in Fig. [2j From the figure one can see that for msusY > 200GeV the expansion 
agrees well with the exact result showing a relative deviation below 10%. For all ttisusy > 
m t the relative correction to the threshold cross section stays below 0.6%. The size of the 
SQCD corrections in a non-universal SUSY mass scenario is shown in Fig. [3]where A S( ^ CD 
is plotted as a function of m^ 2 and rrig for = m t . Since our results are 7r-periodic 
in 6i we have chosen for illustration the four values % G {0, j , ~, ^}. The figures show, 
that only for light masses of the second top squark (mj 2 < 2m t ) A S( ^ CD can have a strong 
dependence on m^ 2 . 

In general one observes corrections below 1% which become negligible for large masses 
of the SUSY particles. A correction factor above 1% is only observed for 9$ = 7r/4 and 
relatively light gluino masses of the order of the top quark mass which are excluded within 
the MSSM [30]. 

4 Electroweak corrections in the THDM and the 
MSSM 

QCD corrections only affect the jit/Zti vertex. On the other hand, electroweak cor- 
rections require also the inclusion of the e + e~ , y/e + e~Z vertex and furthermore of gauge 
boson self energy and box contributions which are necessary in order to arrive at a finite 
and gauge parameter independent result. Typical Feynman diagrams contributing to the 
individual building blocks are shown in Fig. H] for the SM and in Fig. [5] for the MSSM. 
Due to the renormalization procedure (we follow Ref. [31,32]) also W boson and fermion 
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Figure 2: A S( ^ CD = 2dg computed from Eq. ( fT6l) as a function of ttt-susy- The solid line 
represents the exact result and the dashed curve the expansion including terms up to 
order {m 2 t /ml VSY ) 5 . 

selfenergy contributions have to be computed which are also shown in Figs. H] and [51 
They are used in order to render the four building blocks individually finite which is quite 
convenient to deal with. 

In a first step we have used our set-up in order to compute the SM contribution. We 
find complete agreement with Refs. [18,19,21]. Afterwards the THDM model has been 
considered and the results from Ref. [20] have been reproduced^ Let us note that for 
vanishing bottom quark mass the corrections in the THDM model can easily be obtained 
from the analytical results for Higgs- and Goldstone boson contribution calculated in the 
SM by adjusting the coupling factors and boson masses in the loop diagrams. 

Results for the MSSM are not yet available in the literature. However, it is possible to 
compare our results for the vector boson selfenergies with Ref. [33] where the top quark 
production has been considered above the threshold. As far as the box contribution is 
concerned, new kind of diagrams occur in the MSSM where the electron and positron in 
the initial state are not part of the same fermion line (and similarly for the top quarks 
in the final state), cf. Fig. [5](a). Due to the different tensor structure, originating from 
the Majorana character of charginos and neutralinos, it is not straightforward to process 
these contributions with our set-up. On the other hand, it is possible to extract the 
relative correction to the cross section at the threshold by taking the limit s — > 4m^ since 
these diagrams only involve heavy particles inside the loop. However, due to the numerical 
properties of the loop functions [24] the limit can not be taken naively. Instead we evaluate 
the result of Ref. [33] for the box contribution above threshold and extrapolate to s — 4mf. 

2 In Ref. [20] the expression for af^o is proportional to cos(/3— a) which should be replaced by sin(/3— a). 
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(c) A SQCD (6~ t = 7T /2) (d) A SQCD (0 £ = 3tt /4) 

Figure 3: A S( ^ CD as a function of m~ t2 and m- g (normalized to the top quark mass) for 
m £l = rrit and different values of the mixing angle 8^. 

In this way we obtain the threshold contribution for the new box diagrams with three 
significant digits which is sufficient for the phenomenological analysis. We have applied 
the same procedure for the SM box contributions which provides both a cross check on our 
analytical calculation and the very procedure for extracting the threshold contribution. 

Due to the occurrence of many different masses and mixing angles the remaining 
general expression is quite lengthy in the case of the MSSM. Thus, in the following we 
will only discuss the numerical effects. In Ref. [34] a package is provided which allows the 
numerical evaluation of the corrections described in this paper. It uses Mathematica as 
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Figure 4: Typical Feynman diagrams contributing to A 



SM 




Figure 5: Typical Feynman diagrams contributing to A 



MSSM 



front-end and calls Fortran for the time-consuming parts of the calculation. In addition 
an interface to SPheno [35] is provided, which generates numerical values for the masses 
and mixing angles on the basis of a certain SUSY breaking scenario. 

In the numerical discussion we will restrict ourselves to the SUSY breaking scenario 
based on minimal supergravity (mSUGRA) and use the Snowmass Points and Slopes 
(SPS) [36, 37] in order get an impression of size of the corrections. In addition to the 
five mSUGRA parameters mo, m^, tan/3, Aq and sgn(/i) (cf. Tab.l) which serve as 
input for the spectrum generator we use the following input values for the remaining SM 
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Points 


Slopes 


Label 


m 


mi/2 


A 


tan (3 


m 


A 


SPSla' 


70 


250 


-300 


10 






SPSla 


100 


250 


-100 


10 


0,4m 1/2 


-0, 4m 1/2 


SPSlb 


200 


400 





30 






SPS2 


1450 


300 





10 


2m 1/2 + 850 





SPS3 


90 


400 





10 


0,25m 1/2 - 10 





SPS4 


400 


300 





50 






SPS5 


150 


300 


-1000 


5 







Table 1: Input values for the SPS scenarios as defined in references [36,37]. All masses 
are given in GeV and sgn(/i) = 1. 





SPSla 


SPSla' 


SPSlb 


SPS2 


SPS3 


SPS4 


SPS5 


^SM EW 


0.152 


0.151 


0.149 


0.149 


0.149 


0.150 


0.149 


^THDM EW 


0.097 


0.096 


0.093 


0.091 


0.093 


0.099 


0.094 


A MSSM EW 


0.096 


0.096 


0.093 


0.089 


0.093 


0.101 


0.094 



Table 2: Numerical values for A XEW l6 {SM, THDM, MSSM} for variouse SPS scenar- 
ios. 

parameters [38-4oJf| 

m w = 80.40 GeV, m z = 91.1876 GeV, c 2 w = m 2 w /m 2 z , 

m t = 173.1 GeV, m b = 4.2 GeV, a' 1 = 137.036, 

Aa hld( m z) = 277.45 x 10~ 4 , A«i ep (m z ) = 314.97 x 10~ 4 . (17) 

In a first step our Mathematica program transfers the input values to the spectrum 
generator SPheno [35] which produces numerical values for all unknown MSSM parameters 
relevant for our analysis. The output is automatically imported into Mathematica and 
afterwards used in order to evaluate the THDM or MSSM corrections. More details 
about the functionality of our package is provided via the usual Mathematica internal 
documentation and example files which in addition automatically generate the plots and 
tables shown in this paper. 

The numerical impact of the corrections in different mSUGRA scenarios can be seen 
in Tab. Swhere A SMEW , a thdmew and A MSSMEW are eva i uate d f or seV eral SPS pointfl 
Note that A EW varies since the SM Higgs boson is identified with the lightest MSSM 
Higgs boson. 

3 Following Ref. [41] we replace light fermion contributions to the derivative of the photon vacuum 
polarization function by Aa^ d (mz) and Aai cp (mz)- 

4 We add EW to the superscript in order to make clear that only electroweak and no strong corrections 
are considered 
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The SM corrections amount to a sizeable shift of about 15% which get reduced by 
roughly 5% to 6% in the case of the THDM. The main reason for this reduction is the 
smaller coupling of the top quark to the light Higgs boson. At the same time only 
numerically small contributions arise from the diagrams involving heavy Higgs bosons. 
In Tab. [2] one observes only a marginal difference between the THDM and the MSSM. 
It is thus instructive to have a closer look at the depedence on mi/ 2 as suggested by the 
SPS scenarios. For illustration we show in Fig. [6] the comparison of A SM EW , A THDM EW 
and A MSSM EW for SPS1 and SPS2. In both cases we observe only small corrections 
beyond the THDM, i.e. from the neutralino and chargino sector of the MSSM. Larger 
deviations of the order of 0.5% are only observered for those values of 7711/2 where the 
corresponding chargino masses are close to the top quark mass. This becomes clear in 
Fig. 7 where we show the correction of the finite building blocks separately for the case 
of SPSla. One can see a relatively strong variation in the dashed curve which shows 
the contributions from the charginos. It is interesting to note that the peak around 
777.1/2 ~ 150GeV in Ab ox is clearly visible in Fig. O (a) whereas a cancelation among the 
various parts occures for the peak close to mi/2 ~ 250GeV. For comparison we plot in 
FigfT^d) the contribution from SQCD (dashed-dotted). For m 1( / 2 > 150GeV it is smaller 
than the chargino contribution. Corrections above 0.5% are only reached for relatively 
small values of mi/ 2 which corresponds to small values of the gluino mass. 
Let us finally mention that we performed our calculation in the THDM and MSSM for 
finite bottom quark mass and investigated possible large corrections for higher values of 
tan j3. However, even for the mSUGRA scenario SPS4 where tan f3 = 50, the result for 
massless bottom quark is Ajjr^f EW = 0.099, thus the effect of finite bottom quark mass 
adds 0.002 to A^ S £M ew ( see Tab. M>. 



5 Conclusions 

In this paper we investigated the complete weak and strong one-loop corrections within the 
MSSM to the threshold production of top-quark pairs at a future e + e~ linear collider. For 
the SM, THDM, QCD and SQCD corrections we confirmed the results in the literature, 
the genuine supersymmetric electroweak corrections are new. 

As far as the numerical importance is concerned, the electroweak SM corrections 
amount uo to +15% for light Higgs masses. After extending the Higgs sector we ob- 
serve for the SPS scenarios a screening of about —5% to —6% in the THDM (type II). 
The pure supersymmetric corrections from the chargino, neutralino and the strong sector 
are below 1% in most of the parameter space. 
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A On-shell counterterms 



In this appendix we discuss the definition of the counterterms appearing in Eq. ( [131) . In 
order to define the fermionic on-shell counterterms one needs the coefficient functions of 
the tensor decomposition from the corresponding fermion selfenergy: 

E(q, m) =m E s (q 2 , m) + i E v (q 2 , m) + j^ 5 E a (q 2 , m). (18) 

They can be extracted with the help of the following projections: 

-tr{E(q,m)}, 

i 

;tr{$E(q,m)} , 
-tr (7 5 ^S(g,m)} . 



E s (q 2 ,m) 


1 


Am 


E v (q 2 ,m) 


1 

" Aq 2 


E a (q 2 ,m) 


1 



The wave function counterterms are then given by 



5Zl 



6Z{ 



E a (m 2 f ,m f ) . 



E v (q 2 ,m f ) + E s (q 2 ,m f ) 



(19a) 
(19b) 

(19c) 

(20) 
(21) 



B Mixing matrices 

Let us for definiteness provide in this Appendix the definition of the mixing matrix = 
(Uij) used in Eq. [TD and [151 We work with flavor diagonal sfermion mixings, where the 
left and right handed sfermion fields fi and f R are connected to the mass eigenstates f\ 
and /2 via 

The 2x2 mixing matrix Uj diagonalizes the mass matrix of the corresponding sfermion 
/: 

U/m 2 -Ut = ( h ^ ] . (23) 



J2 



In the case where the mass matrix rn 2 - contains only real entries, one can choose to be 
orthogonal. For its parameterization only one angle 9j is needed and the transformation 
from mass to gauge eigenstates can be written as follows: 

fi = h cos 9 f - + f R sin 9 f - , 

h = f R cos9 f - f R sm9 f . (24) 
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